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Background 


Experimental techniques for measuring the kinetic behavior of liquid-liquid droplet 
dispersions held isothermally in thin cuvettes were developed at the Space Science 
Laboratory, Low Gravity Research Division, NASA Marshall Space Flight Center. The 
observational method, based on optical holography, was developed originally by Dr. D.O. 
Frazier and his scientific staff, including W.K. Witherow, B.R. Facemire, J.P. Downey, 
and J.R. Rogers. This holographic data set at MSFC is the largest set of quantitative 
coarsening data available. The Principal Investigator (M.E. Glicksman) and the Co- 
Principal Investigator (V.F. Fradkov), both at Rensselaer Polytechnic Institute, were asked 
originally by NASA to assist scientists at MSFC in interpreting the observed coarsening 
kinetics as the liquid-liquid system (isopycnic phases consisting of succinonitrile-rich 
drops in a water-rich matrix) was annealed over long periods of time. 

Cooperative research activities originally initiated under NASA contract NAG8-864 
included developing novel numerical techniques for predicting the lifetimes of individual 
droplets located within a specified area of a hologram. These predictions were based on 
solving the multiparticle diffusion problem in two spatial dimensions, subject to certain 
boundary conditions applied to the droplets as 3-dimensional objects. This so-called mixed- 
dimensional interpretation of the coarsening kinetics, as observed by Frazier et al., proved 
extremely beneficial in correlating the holographic observations with a theoretical 
framework. This collaborative research work between MSFC and Rensselaer appeared 
subsequently as the following two archival papers: 

1) J.R. Rogers, J.P. Downey, W.K. Witherow, B.R. Facemire, D.O. 

Frazier, V.E. Fradkov, S.S. Mani, and M.E., Glicksman, “Coarsening of 
Three-Dimensional Droplets by Two-Dimensional Diffusion: Part I 
Experiment,” Journal of Electronic Materials , Vol. 23, No. 10, pp 999- 
1006, (1994). 

2) V.E. Fradkov, S.S. Mani, M.E., Glicksman, J.R. Rogers, J.P. 

Downey, W.K. Witherow, B.R. Facemire, and D.O. Frazier, “Coarsening 
of Three-Dimensional Droplets by Two-Dimensional Diffusion: Part II 
Theory,” Journal of Electronic Materials , Vol. 23, No. 10, pp 1007-1013, 

(1994). 
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Annual Progress 


The Principal Investigator and Co-Principal Investigator have developed a novel 
numerical method for performing multiparticle diffusion calculations that permit statistical 
analysis of the NASA/MSFC droplet holograms. The new method uses instantaneous 
“snap-shots” calculations of many diffusionally interacting particles (up to several 
thousand). Heretofore, time-marching techniques were needed to determine the coarsening 
rates in mutiparticle clusters. Such time-dependent calculations were so computationally 
laborious that the kinetic behavior of only a few hundred particles could be determined even 
with the largest computers available. Now clusters containing thousands of particles can be 
assessed numerically with only modest computing power. 

The theoretical basis for the new method, based on linear perturbation analysis, 
along with some initial results obtained for three-dimensional clusters has been published 
recently. This paper, attached to this report in pre-print form, appears as Appendix 1 , is 
entitled: 


“Coarsening kinetics in finite clusters,” Physical Review E, Vol. 53, No. 

4, pp. 3925-3932, (1996). 

The new method, detailed in Appendix 1, specifically permits assessment of the 
diffusion screening distance (Debye length) in two-phase mixtures, such as being studied 
by Frazier et al. at MSFC. The availability of the MSFC hologram data base for droplet 
coarsening can be analyses further by applying the methods disclosed herein. 

The Co-Principal Investigator, Dr. Fradkov, is currently spending a 12- week 
period of the Summer, 1996, at MSFC, working with the NASA scientists in Dr. Frazier’s 
Branch. Further details of the application of these methods by MSFC scientists are being 
worked out during the visit to Huntsville, Alabama, by Dr. Fradkov. The Principal 
Investigator also plans at least one trip to MSFC during the Summer of 1996 to interact 
directly with the MSFC scientists. 


Future Plans 

If the current cooperative agreement is funded for 1997 and extended for one year 
beyond the requested extension date of September 30, 1996, the following plan of research 
activities will be followed: 
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• Multiparticle calculations will be provided for studying higher volume fractions 
than investigated for MSFC during 1996. Our initial study, as detailed in Appendix 1, was 
limited to a maximum volume fraction (in three dimensions) of 0.01. The Principal 
Investigator feels that this restriction can be relaxed by using novel techniques that extend 
these numerical computations to “dipolar” order, as compared to the “monopolar” 
computations accomplished to date at low volume fractions. This approach is especially 
valid for analyzing the MSFC droplet coarsening hologram data base, which require 
solving the diffusion equation in only two space dimensions. The equation set for solving 
the diffusion field for large two dimensional clusters is not as difficult to invert as is the one 
for a three dimensional cluster. Implementing this plan will require cooperative decisions 
with the cognizant NASA scientists at MSFC. 

• The Principle Investigator would attempt to implement some fast inversion 
methods for solving the “snap shot” diffusion equation to dipolar order. If successfully 
accomplished this development would lead to coarsening kinetic coefficients useful for 
analyzing data obtained at higher area fractions in two dimensions. This element of our 
future work would be carried out cooperatively with applied interested mathematicians here 
at Rensselaer Polytechnic Institute, at the Courant Institute, New York, NY, and at the 
National Institute of Standards and Technology, Gaithersburgh, Maryland. 
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Abstract 

We address the problem of diffusional interactions in a finite sized cluster of spherical 
particles for volume fractions, V v , in the range 0-0.01. We determined the quasi-static 
monopole diffusion solution for n particles distributed at random in a continuous matrix. 
A global mass conservation condition is employed, obviating the need for any external 
boundary condition. The numerical results provide the instantaneous (snapshot) growth 
or shrinkage rate of each particle, precluding the need for extensive time-dependent 
computations. The close connection between these snapshot results and the coarse- 
grained kinetic constants are discussed. A square-root dependence of the deviations of 
the rate constants from their zero volume fraction value is found for the higher V v 
investigated. This behavior is consistent with predictions from diffusion Debye-HUckel 
screening theory. By contrast, a cube-root dependence, reported in earlier numerical 
studies, is found for the lower V v investigated. The roll-over region of the volume fraction 
where the two asymptotics merge depends on the number of particles, n, alone. A 
theoretical estimate for the roll-over point predicts that the corresponding V v varies as ri 2 , 
in good agreement with the numerical results. 
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I. Introduction 

Diffusionai coarsening represents an important and commonly observed kinetic 
process in microstructural evolution. Coarsening can occur among several microstructural 
constituents ranging from the primary phases to widely dispersed precipitates. The wide 
range of volume fractions, V v , encountered in phase coarsening, makes it essential to have 
a fundamental theory of coarsening which treats the volume fraction as an input parameter 
of the microstructure. The classic coarsening theory developed by Todes [1] and Lifshitz 
and Slyozov [2] (TLS theory) is limited to infinitesimally small volume fractions. To date, 
scores of papers have been published attempting to extend TLS theory to finite volume 
fractions intrinsic to real microstructures. Several excellent reviews have been published 
surveying these models along with relevant experiments [3-5]. Nonetheless, even finding 
the initial corrections to TLS theory in the limit of small volume fractions remains an open 
question. Specifically, the analytical theory for infinite systems suggested by Marqusee 
and Ross [6] and Marqusee [7], based on the diffusion analog of the Debye-Htickel 
screening effect, leads to deviations from TLS in the coarsening rate, which are 
proportional to V v . By contrast, several numerical analyses [8,9] dealing with finite 
systems clearly show that this deviation is, instead, proportional to V v m . The purpose of 
the present paper is to analyze how the coarsening rates change in a finite cluster of a 
specified particle density. We will show numerically that both V v m and V v m behaviors 
may occur, depending on the number of particles comprising the cluster and its volume 
fraction. The crossover between these behaviors will be derived analytically using the 
Debye-Hiickel screening concept. 

II. Theoretical Background: Mean-Field Approaches 

The simplest model of coarsening in infinitesimally sparse systems [1] assumes that 
all precipitate particles remain spherical and remote from each other, and that their positions 
remain stationary. The description of isothermal coarsening is given in [1] in terms of the 
radius distribution function, F(R,t), where R and t are the particle radius and the time. The 
norm of F(R,t) is based on the total number of particles per unit volume, N v , that is 

O* 

fF(*,/)d* = iV (1) 

0 

With this normalization, the volume fraction, V v , of the dispersed phase is given by 

V.^FiRj^R'dR^^-N^). ( 2 ) 

o - 1 

Supersaturation of the matrix solution is assumed to be small, so that nucleation of 
new particles is precluded. The continuity equation for particles in size space is 

(3) 
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where v(/?) is the time-rate of change of the radius, R, of a particle. To obtain v(R), 
further simplifying assumptions were made: (a) the kinetics of coarsening is limited by 
volume diffusion; and (b) the diffusion is quasi-static. Therefore, the diffusion equation 
describing the concentration field, c(r), in the matrix phase reduces to the Laplace equation 


V 2 9 = 0, (4) 

where <p = (c-c 0 )/c 0 and c 0 is the equilibrium solubility. The boundary conditions at the 
surface of the i-th particle are specified through the Gibbs-Thomson local equilibrium 
relation 


X 


( 5 ) 


where X is the capillary length given by 


X = 


_2yQ 


k B T 


( 6 ) 


Here, y is the particle-matrix specific interfacial energy, Q is the atomic volume of the 
dispersed phase, k B is the Boltzmann constant, and T is temperature. Additionally, 
anisotropy of the interfacial energy and strain effects were ignored. 

The solution to eq. (4), subject to the boundary conditions (5) at the interface of 
each particle, may be written in the form of the Coulomb potential 1 , 

XB i 






+ <P-> 


( 7 ) 


where <p_ is the background matrix diffusion potential, assumed uniform throughout the 
matrix phase; r is a field point always in the matrix; r, locates the center of the i-th particle 
The diffusion analog of electrostatic charge is the volume flux per steradian, B ; , given by 

R 
R' 


£,= 


where the critical radius, R\ is given by 


D ^ 

1_*L 

*) 


( 8 ) 


R'=- 


<P« 


( 9 ) 


The time-rate of change of the radius, v(R), used in eq. (3), is directly connected to B i as 




2 XDc^l 


*1 


B, 


( 10 ) 


Equation (8) provides the values of B p but the value of (p^ remains unknown. Todes 
suggested using the global conservation of volume of the dispersed phase, neglecting 
changes in the supersaturation within the matrix phase [1]: 


1 The solution given by eq. (7) is a monopole approximation, where the particles are treated as point 

sources and sinks. This is justified by the large distances between the particles compared to their radii. 
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( 11 ) 


o* «• 

— (V( /?,/)• — R 3 dRs [^-F(R,t)-—R 3 dR=0. 
dt) K, 3 J dt K ’3 

0 0 

Using the continuity equation to substitute for dF/dt, one obtains 

om 

\F(R,t)dR 


<P-=:f 


a m- orR ' = ( R ) 

\F(R,t)RdR ' ' 


( 12 ) 


(14) 


(15) 


Todes also introduced a dimensionless variable p =R/R‘. The dimensionality of the 
diffusion-limited coarsening process implies that 

R = oc(2 ‘kDc^Qt} (13) 

where a is a dimensionless kinetic coefficient to be defined later. Todes then found a one- 
parameter family of self-similar solutions for the distribution function, F(R,t), given by the 
product-function form 

With eq.(14), the continuity equation, (3), reduces to the following O.D.E. 

df ~4p 3 +(6/a 3 )(l-p) dp 

/ (3/<x 3 )(l-p) + p 3 P ' 

The number of real roots of the denominator, <D=(3/a 3 )(l-p)+p 3 , on the r.h.s. of eq.(15) 
depends on the value of a, which in turn influences the qualitative behavior of the solution 
toeq.(15). Figures la-c show the three possible behaviors of the denominator. For 
a>{213) 2ri , figure 1-a, there is no real root of the denominator, and the solution to eq.(15) 
is analytic up to infinity. This solution, however, causes a logarithmic divergency of the 
integral for the total volume of the particles. Truncating the logarithmic tail, however, 
destroys self-similarity of the distribution function, ftp), because the relative density of 
particles at large p values is gradually exhausted. For a<(2/3) 2/3 , figure 1-c, the solution 
must be restricted between zero and the first positive root of the denominator and set equal 
to zero beyond, in order to prevent singular behavior of the distribution function at the 
roots of the denominator. These solutions are, however, unstable with respect to 
appearance in the system of a particle outside of the permitted range as a result of 
coalescence. Lifshitz and Slyozov [2], therefore, stated that the only stable solution to the 
continuity equation, (3), corresponds to such a value of a that makes the two positive roots 
of the denominator coincide, figure 1-b. That implies 

a 7-tS = ( 2/3 ) t (16) 

It is important to note that although the TLS solution to coarsening kinetics does not 
depend formally on the value of the volume fraction, V v , it is not, in fact, extensible to any 
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non-zero volume fraction. The restriction of the TLS solution to zero volume fraction is a 
consequence of the contradiction between the infinite extent of the Laplace diffusion field 
and the total neglect of direct interactions among the particles. This contradiction manifests 
itself even stronger in two dimensions, where an analogous TLS theory cannot even be 
developed without imposing some arbitrary cutoff radius to the diffusion field, due to its 
logarithmic far-field behavior for each individual particle [8, 9]. 

There are two general approaches used to incorporate the influence of volume 
fraction on coarsening kinetics. The earliest approaches were based on arbitrary restriction 
of the extent of the interparticle diffusional interactions to the average interparticle 
separation (see [17] for a review). In sparse systems the nearest neighbors, of course, 
cannot shield the Laplacian field. Theories based on models using interparticle distance for 
the cutoff are generally in poor agreement with experiment and numerical simulations [3, 4, 
17]. For higher volume fractions, 0.1<V\,<0.6, however, direct screening by nearest 
neighbors becomes more reasonable, but generally the models that use a single cutoff 
radius do not provide meaningful results for the specified range of volume fractions (see 
ref. [10] for a review). 

Marqusee and Ross [6] and Marqusee [7], by contrast, restricted the extent of the 
Laplacian diffusion field by taking into account screening effects in active two-phase media 
containing a distribution of diffusional sources and sinks. Instead of the Laplace equation, 
Marqusee and Ross obtained the Poisson equation for the spatially coarse-grained 
backgound diffusion potential, <p, namely 

V 2 (p = -4rto, (17) 

in which the source/sink density, a, in the matrix space surrounding a particle is given by 

a(r) = JX5 • F(R,t) d R = J(X - R<p(r))F(R,t)d R. (18) 

0 o 

Note that here the background potential in Eq.(8), (p Ml is replaced by the coarse-grained 
potential, cp, consistent with Debye’s approach. Carrying out the indicated integration on 
the right-hand side of eq. (18) gives 

a(x) = N V X - N v ( J?)< p(4 (19) 

where N v and ^R^ are defined by eqs. (1) and (12), respectively. Substituting eq. (19) 
into eq. (17) yields the diffusion analog of the Debye-Hiickel equation, namely 



V 2 (p-K 2 (cp — <P-) = 0, 

(20) 

where 


1^=4 1tNy(R), 

(21) 

and 


II 

8 

9- 

(22) 
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Equation (20) is well known from the theories of ionic plasmas and electrolytes, 
where it is used to describe the electrostatic field developed in an active medium containing 
a distribution of quantized charges. In these theories, the second term on the l.h.s. of eq. 
(20) describes the space charge density caused by spatial redistribution of the charges, in 
order to minimize the electrostatic energy. In the case of diffusional coarsening, however, 
the positions of the “charges” (diffusional source/sink strengths) are fixed, but their values, 
the 5,’s, depend on the local background potential, cp(r), consistent with eq. (8). 

The solution to eq. (20) subject to the boundary conditions, eq. (5), is in the form 
of a Yukawa potential 

^( r ) = Sr~i ex p(" K: 'l r_:, 5'l) +( P-’ (23) 

1 r ^1 

where now the B' s are defined as 


5 ,= 




(24) 


denotes the dimensionless volume fluxes (per steradian) for particles in an 
infinitesimally dilute system, as given in the TLS theory by eq. (8). The coefficient k is the 
reciprocal of the diffusional analog of the Debye screening length. It represents the natural 
cut-off distance for direct particle interactions via the diffusion field, beyond which the 
particles are isolated from each other by the intervening two-phase medium. As seen from 
eq. (24), k /?, expressly captures the deviations in the rate of particle size evolution from the 
TLS limit. It is possible to use the same global stability analysis introduced by Lifshitz and 
Slyozov [2] to obtain the characteristic kinetic constant, a, which determines the overall 
rate of coarsening as defined in eq. (13) (see [6] for details): 


a(F K ) = a m (l + 0.740 f; /2 + H.O.T.), (25) 

where a ns =(2/3) 2n is the TLS-limit for infinitesimally diluted systems (V v =0). Equation 
(25) shows that the average particle size grows more rapidly than in the TLS limit. 
Moreover, this correction, reflecting diffusional Debye screening, implies other first-order 
corrections, including the broadening of the size distribution function and the difference 
between the average and the critical radii; all these corrections are proportional to the 
square root of the volume fraction. 


III. Kinetic Properties of Finite Coarsening Systems 

There are four physically distinct metric scales in any finite coarsening cluster, i.e., 
a finite number of particles confined to a finite volume: 1) the size of the particles, as 
characterized by A R ~ (ft)-, 2) the interparticle spacing, as characterized by A v = Nf n , 
where N v is the total number of particles per unit volume, eq. (1); 3) the Debye screening 
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radius, A 0 s k~‘; and 4) the extent of the total coarsening system. A,,*, as defined for a 
“spherical cluster” by A tot = (4 tc/3) * 3 (n/ iV,,)' 73 , where n is the total number of particles 
remaining at a given time. Note, it is always assumed that the the cluster is locally 
homogeneous and spatially isotropic, and that the overall size of the cluster. A,*, is time 
independent. These assumptions permit one to establish an analog of the volume fraction 
for a finite cluster of particles. 


Vr- 




( 26 ) 


where the denominator in eq. (26) provides an estimate of the spherical cluster volume . 

By contrast with infinite systems, finite systems lack uniform asymptotics, because 
any finite system ultimately achieves its end state as a single particle in finite time. Finite 
systems, instead, may only exhibit intermediate asymptotics. By this we mean that the 
coarse-grained average over time of these cluster parameters becomes identical to those in 
self-similar infinite systems, even though the instantaneous values of the time-dependent 
parameters of finite clusters may fluctuate. All continuum coarsening theories, however, 
exclusively yield quantities that are equivalent to the coarse-grained averages over time for 
finite clusters. Specifically, therefore, the coarse-grained intensive length scale parameters, 
viz. A^ A r , and A D , are expected to scale with time as predicted by continuum coarsening 
theories. 

Weins and Cahn [11] were the first to estimate the diffusional interactions among a 
few particles comprising a small coarsening cluster. Later Voorhees and Glicksman [12, 
13] extended Weins and Cahn’s analysis to hundreds of particles by using periodic 
boundary conditions, and adapted Ewald sums to improve the convergence of the 
numerical calculations. In finite systems, coarsening always occurs via two distinct but 
related processes: the continuous evolution of the particle radii, and the discontinuous 
(instantaneous) disappearances of the smallest particles. The coarse-grained time averages 
depend on both processes. Specifically, during the continuous process, the average size of 
the particles decreases (along with the total interfacial area), whereas the instantaneous 
disappearances cause the average size to increase. The actual evolution of the average 
radius consists of smoothly falling segments connected by vertical jumps whenever a 
particle disappears, as observed in computer simulations [13, 15]. The coarse-grained 
average changed with time as r 1/3 , according to the TLS prediction. 

A recent simulation and experimental study of coarsening in a finite clusters of 
three-dimensional particles interacting through two-dimensional diffusion field [14,15] 
shows similar discontinuous behavior of the average radius, whereas the coarse-graining 
behavior of the average particle radius exhibited a power law of t m , characteristic of the 
continuum coarsening prediction for the mixed-dimensional case (see figure 2). Thus, 
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experiments confirm and simulations show that the coarse-grained intermediate asymptotic 
behavior of real, i.e. finite, coarsening systems is consistent with theoretically predicted 
asymptotic power laws from continuum theories. 

Normally, the numerical data required for coarse-grained averaging are obtained by 
time-dependent computer simulations. Following such procedures for large clusters (tens 
of thousands of particles) is computationally intensive. There would be a great advantage 
in large finite systems to limiting the numerical simulation to a “snapshot” analysis of the 
rate of change of <R> along the continuous segments. It is unclear, however, how the 
instantaneous slopes, v(/? f ), defined by eq. (10), are related to the coarse-grained average 
time dependence of <R>. One can show, using first-order perturbation theory, that the 
relative change in the kinetic constant, a, from the Lifshitz-Slyozov stability analysis is 
proportional to the relative change in the growth rate, v(/? mai ), for the largest particle within 
the TLS -distribution, induced by any small perturbation. Specifically, if a small 
disturbance to the particle growth rates, £(/?), is introduced as 

v(i) = v m («)[l + e(*)], (27) 

where Vju is the non-perturbed value given by the TLS theory, eq.(10). The denominator, 
©, in eq. (15), becomes 

® = p 3 - P(p-l)(l+e). (28) 

where p = 3/a 3 . Following the procedure set out in the earlier section. Theoretical 
Background , we find the eigenvalue for p and the maximum value of p by solving the 
equation set <D=0 and <D=0. Seeking first-order corrections to the eigenvalue, we write 


P — Ptls + ^P 

(29) 

and 


P = Ptls + $P- 

(30) 

where Ptls = 3/2, and P TLS =27/4. 
Keeping only first-order terms, we find 


<D=-5p/2- (27/8)6 = 0 

(31) 

where £ is evaluated at p ^=3/2. This results in 8P = - (27/4)e, and 


a = a 7is l + (3/l6)' /3 £(i^)j. 

(32) 


Here, =(3/2 )<R> and a TLy =(2/3) 2/3 are the non-perturbed values given by the TLS 
theory. For example, if the perturbation, £(/?), is caused by diffusional Debye screening, 
eq. (24), then eq. (32) transforms to eq. (25) from the analysis by Marqusee and Ross 
[7]. 

The second equation, <1 1 =0, may be written as 3<ZV38p=0, or, keeping first-order terms, 

Tf — 6pjLs^P ” SP - Ptl S £ ■ Ptls(Ptls“ 1 = 0. (33) 

This gives the correction in p^, which however depends not on £ but on £'=d£/38p. 
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5p = (3/8) e'. (34) 

The snapshot numerical analysis, to be discussed below, is capable of providing the 
e(R) values for each particle in the cluster. In a finite cluster, however, one does not have a 
statistically valid representation of enough particles at to guarantee an accurate estimate 
of 6(1?^). Consequently, we will average the perturbations over the entire system of 
particles, using the assumption that the global average, <£(/?)>, reflects the behavior of 
e(R m „)• TW S assumption is a rather likely outcome of the self-similar behavior we are 
discussing. 

IV. “Snapshot” Numerical Analysis 

The mathematical formulation for the present numerical study is based on the 
monopole approximation to the quasi-static diffusion solution shown in eq. (7). The 
particle growth rates were calculated for each individual particle in the cluster according to 
eq. (7), using the exact monopole expressions for the dimensionless diffusional 
source/sink strengths, B's, namely 



Here n is the total number of particles in the cluster, A. is a capillary length defined in eq. 

(6), and the expression in the parenthesis represents the local analog of the background 
diffusion potential used in the TLS theory, eq. (8). Equation (35) represents a set of n 
linear equations with n+l unknowns, viz., the B’s and (p... One additional equation is 
needed to close the system. In previous simulations such an equation was provided by an 
external boundary condition, e.g., the periodicity condition used in [9, 13]. In the present 
work we employ, instead, a conservation equation for the total volume of the particles, as 
suggested in [11, 15]: 

= 0. (36) 

A version of the Gauss-Seidel iterative algorithm was employed for solving the set of 
equations given by eqs. (35) and (36), which equations may be written in matricial form as 
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( 37 ) 


A linear projection to nxn subspace was employed to eliminate the last row and column and 
make the matrix, eq. (37), suitable for applying the classic Gauss-Seidel linear solver. 

The cluster parameters in the “snapshot” simulation are: the set of n particle center 
coordinates, { r, } , and the set of n particle radii, {/?,}, both of which are measured in units 
of X. The n particle center coordinates in this simulation are chosen randomly inside a 
sphere of radius A^. Each particle was assigned a radius randomly distributed around the 
average value <R>=A r according to various distributions. Specifically, we used the TLS 
particle size distribution and a rectangular distributions of various widths. An influence of 
the particle size distribution on the kinetic results was not observed. The particles are not 
allowed to overlap. This restriction results in some correlation between particle sizes and 
positions, which tends to vanish with decreasing volume fraction. The volume fraction 
itself was varied by changing the ratio A 3 R / A 3 tot , according to eq. (26), while maintaining 
fixed relative positions and radii of the particles. Specifically, we proportionally changed 
all the interparticle distances and, consequently, changed the size of the cluster and the 
volume fraction, whereas the particle radii remained constant. 


V. Results and Discussion 

Figure 3 shows the dependence on the volume fraction of the observed relative 
deviations, 

Q)) (3f 

( v ( s '°)> 

in the values of the rate of change of the average radius, (v^R,V tr 'jj, compared to that for 

zero volume fraction, (v(R, 0 ]^ . The values at non-zero volume fractions were obtained 

from the snapshot calculations, using eq.(37), whereas the zero volume fraction term is 
obtained from the TLS expression, eq. (8). The results are shown for clusters containing 
10, 100, and 1000 particles. The deviations are plotted on log-log coordinates vs. the 
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volume fraction, V v , given by eq. (26), which was varied for each calculation by changing 
the ratio / A 3 tot , as described above. There appear to be two distinguished limits for the 
slopes, representing the exponents of the leading-order corrections introduced by the 
volume fraction. For the higher ranges of volume fraction the exponent is 1/2, in 
agreement with Debye screening theory, eq.(24). Combining eqs. (2), (8), (10), (21), and 
(24) with the definition of A, eq.(38), we obtain the expression for A in terms of the 
volume fraction. 


A = 


1 -H, 


v 3 V ' 2 




v;“ = 0.695 Vy 


1/2 


3 J 


where the fi,’s are the ith moments of the TLS distribution defined as 


( 39 ) 


OO 

j>F(5,f)d5 

(r)‘n v 


( 40 ) 


The coefficient on the r.h.s. of eq.(39) is slightly different from that appearing in eq.(25), 
because these coefficients reflect different combinations of the moments of the distribution 
function. 

For the lower ranges of volume fraction the exponent appearing in Fig. 3 is 1/3, 
consistent with earlier simulations [13] and theories employing a cut-off radius based on 
the interparticle separation [16]. The rollover point between the two exponents depends on 
the number of particles, n, in the cluster, so that for clusters with larger n the transition 
occurs at progressively lower volume fractions. To explain the observed rollover behavior 
we substitute eq. (2) into eq. (21). This yields the Debye screening length A D in terms of 
volume fraction V v as 

A D =3-' n ((R > )/(R)) m (K)'' n . (41) 


The spatial extent of the cluster, A. ot , can also be estimated in terms of V v as 

A,„ = (4*/3 f’in/K)'' 1 = n'»(R>)"\v y y m . (42) 

Comparison of eqs. (41) and (42) shows that the Debye screening distance grows faster 
with decreasing volume fraction than does the cluster size, holding other system parameters 
constant. This implies that at some small volume fraction, V'., the Debye screening 
distance, A D , exceeds the cluster size, and then diffusional screening is precluded. 

When screening is no longer possible because of the sparsity of the cluster, then the 
deviations from TLS theory no longer have to be proportional to Vj /2 . Indeed, comparing 
the TLS result, eq. (8), with the monopole expression, eq. (35), one obtains for the 
deviation from TLS: 


55 , = 


A 

x 


» 5+55 

Y ; 


( 43 ) 
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For a sparse system, 5 5 y on the r.h.s. of eq. (43) may be neglected. The r.h.s. of eq. (43) 
may be then estimated, using the fundamental metric scales of the system, as In 

terms of the volume fraction, one may express eq. (35) as 


Ml. a i W - W r> 

B i A" N- m l R >\ ,n ' . (44) 


Equation (44) shows that for extremely sparse clusters, the leading order deviation from 
TLS kinetics goes as V', n . This result is consistent with the asymptotic analysis performed 
in [12] for a periodic system. 

The rollover point at which the exponent changes from 1/2 to 1/3 occurs when the 
Debye screening distance becomes comparable to the cluster size. The corresponding 
volume fraction, Vy , may be estimated from the condition A tot =A 0 , which with the use of 
eqs. (41) and (42) gives the result 



1 

27« 2 " 


(45) 


Figure 3 also shows the theoretical rollover points calculated from eq.(45) denoted 
by vertical arrows for clusters containing n=10, 100, and 1000 particles. The predictions 
from eq.(45) appear satisfactory for n=100 and 1000, despite the fact that the clusters 
investigated here have relatively few particles to provide the Debye screening behavior. We 
note, that a cluster with 1000 particles contains only about 6 interparticle distances from its 
center to its periphery. The cluster containing 10 particles fails to show screening behavior 
for volume fractions approaching 1%. Diffusional interactions in clusters with volume 
fractions above about 1% would not be represented accurately by the monopole 
approximation. 


VI. Conclusions 

1) Monopole snapshot calculations provide an efficient method for estimating the 
coarse-grained rate constants in finite clusters of coarsening particles. Prior to development 
of this technique, fully time-dependent simulations were required to obtain this 
information. 

2) Extending the definition of volume fraction to a finite cluster permits 
investigating the influence of non-zero volume fraction on the coarsening kinetics. 

3) The present results confirm that in finite clusters in the limit of low volume 
fraction the leading correction to the TLS rate constant is proportional to Vy 1 . For higher 
volume fractions, or for a sufficient number of particles, Debye screening occurs and the 
leading-order term becomes proportional to Vy 2 . 

4) The rollover point between the two kinetic behaviors, that is, the critical volume 
fraction at which the Debye screening distance is comparable to the cluster size is given by 
eq.(45) and depends on the number of particles in the cluster. For clusters with only one 
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thousand particles, the rollover point is already as low as ca. 3- 10' 8 . This implies that in 
most experimental coarsening systems Vy 3 behavior should not occur. 

5) Inasmuch as Debye screening is not direcdy applicable to small clusters, its 
apparent success in explaining the results of these calculations require further analysis. 
Larger-scale calculations with clusters containing much greater numbers of particles will be 
required to resolve this interesting issue, as well as to provide additional information on 
correlation effects and the influence of the cluster periphery on coarsening kinetics. 

6) At volume fractions above ca. 1%, particle interactions become more complex 
than can be described with monopole approximations and Debye screening. Other 
techniques that incorporate spatial correlations among the particles, shape changes, direct 
screening, and higher-order multipole interactions must be employed. 
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Figure Captions 

1. Behavior of the denominator, <D, of eq.(15): 

(a) tx>(2/3) 2/3 ; 

(b) a=(2/3) 2/3 ; 

(c) oc<(2/3) 2/3 . 

2. Instantaneous and coarse-grained time dependence of the average particle size observed 
in computer simulation of mixed-dimensional cluster coarsening [14, 15]. 

3. Relative deviation. A, of the kinetic rate constant for coarsening of finite clusters at 
various volume fractions from that at zero volume fraction. These data are obtained by 
snap-shot monopole calculations. Theory refers to eq.(39). Arrows correspond to the roll- 
over volume fractions where the exponent changes from 1/3 to 1/2, as predicted by 
eq.(45). 
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